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The Waxman-Peck theory of the population genetics is discussed in regard of soil bacteria. Each 
bacterium is understood as a carrier of a phenotypic parameter p. The central aim is the calculation 
of the probability density with respect to p, <l>{p,t;po), of the carriers living at time t > 0, provided 
that initially, at to = 0, all bacteria carried the phenotypic parameter po = 0. The theory involves 
two small parameters: the mutation probability jj, and a parameter 7 involved in a function w{p) 
defining the fitness of the bacteria to survive the generation time t and give birth to offspring. The 
mutation from a state p to a state q is defined by a Gaussian with a dispersion am- The author 
focuses attention on a function (y3(p, t) which determines uniquely the function <?(p, i;po) and satisfies 
a linear equation (Waxman's equation). The Green function of this equation is mathematically 
identical with the one-particle Bloch density matrix where fi characterizes the order of magnitude of 
the potential energy. (In the a;-representation, the potential energy is proportional to the inverted 
Gaussian with the dispersion a^.) The author solves Waxman's equation in the standard style of a 
perturbation theory and discusses how the solution depends on the choice of the fitness function 
In a sense, the function c{p) = l — w{p)/w{0) is analogous to the dispersion function E(p) of fictitious 
quasiparticles. In contrast to Waxman's approximation where c{p) was taken as a quadratic function, 
c(p) ~ 7p^, the author exemplifies the problem with another function, c{p) = 7(1 — exp(— ap^)], 
where 7 is small but a may be large. The author shows that the use of this function in the theory 
of the population genetics is the same as the use of a non-parabolic dispersion law E = E{p) in the 
density-matrix theory. With a general function c(p), the distribution function <?(p, t; 0) is composed 
of a delta-function component, N{t)5(p), and a blurred component. When discussing the limiting 
transition for t — > 00, the author shows that his function c(p) implies that N(t) N{oo) 7^ in 
contrast with the asymptotics N{t) ~+ resulting from the use of Waxman's function c(p) ~ p^. 

PACS numbers: 02.50.Ey, 05.10.Gg, 87.10.+C, 87.23. Cc 



I. INTRODUCTION 



During the past century, the development of the quan- 
tum theory has been paralleled with the development 
of the theory of stochastic processes. When comparing 
the time-dependent Schrodinger equation of the quan- 
tum mechanics with the Fokker-Planck equation of the 
stochastic dynamics, we may assert that both these equa- 
tions are of the same kind. Mathematically speaking, 
both these equations are linear second-order partial dif- 
ferential equations of the parabolic type. The Fokker- 
Planck equation can be derived for any process which 
can be described by the Langevin equation ridx{t)/dt — 
F{x{t)) = f{t) (cf. e.g. [1, 2]). In the Langevin equa- 
tion, 77 > is a deterministic constant, F{x) is a de- 
terministically defined real- valued function and f{t) is a 
stochastically defined zero-centered stationary Gaussian 
white-noise function. (As a rule, x{t) and f{t) are con- 
sidered as real random functions of the real time variable 
t. The stationarity of fit) is meant in the stochastic 
sense.) Using the angular brackets ( ) for the averaging 
with respect to the randomness of /(i), we assume that 
(/(t)) = and {f{ti)f{t2)) = A 5{ti - tj) at all time 
instants t, ti, t2. The initial value of x(t), x{0) = xq, 
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is usually a deterministic value given in advance. The 
Fokker-Planck equation concerns the conditional proba- 
bility density P{x,t\xo) = {Six{t) - x)) > 0: 



dP{x,t\xo) A d'^P{x,t\xo) 



dt 



2772 



dx"^ 



ld_ 

J] dx 



F{x)P{x,t\xo] 



P{x, -t-0|a;o) = S{x — xq) 



(1) 
(2) 



If x{t) is the instantaneous position of a Brownian par- 
ticle moving along a line, we may speak of the diffusion 
coefficient D (such that 2ri^D = yl > 0) and of the mo- 
bihty 1/77 of the particle. Then F{x) = — dVg(x)/dx is 
a driving force, f{t) is the Langevin stochastic force and 
(|l|) is the governing equation of the Brownian dynamics. 
Equation can be transformed into an equation of the 
Schrodinger type by the substitution 

Pix,t\xo) = exp [V^ixo)-VQ{x)]^ Rix,t\xo). (3) 

The function R(x,t\xo) obeys the equation 
dR{x,t\xo) A d'^R{x,t\xo) 



dt 
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where 



dx'^ 



dx 



V{x)R{x,t\xo), (4) 



d^V^ix) 



da;2 
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[F{x)Y dF{x) 



A 



(5) 



Equation (Q) is known as the real Schrodinger equation. 
If t is replaced by fip, where k^T = 1//3, equation 
with some change of symbols, becomes the Bloch 
equation for the one-particle canonical density matrix of 
boltzons of a constant (effective) mass in the thermo- 
dynamic equilibrium at the temperature T [3]. Thus, 
when comparing equations (|l|) and (^, we can always 
juxtapose the Brownian dynamics and the quantum the- 
ory. (This juxtaposition can also be based on Feynman's 
path-integral theory [4, 5].) The transition from the for- 
malism of the Brownian theory to the formalism of the 
quantum theory is easy since if we define the driving 
force F{x) = — dVg(a;)/dx, we can directly calculate the 
function V{x) according to equation (||). {Vice versa, if 
one tries to find F{x) corresponding to a given function 
V{x), one has to solve equation (^) which is nonlinear. 
Equation (|^) is known as the Riccati equation. Cf. any 
handbook on nonlinear differential equations — e.g. [6]. 
Recently the usefulness of the Riccati equation in solving 
various problems of classical and quantum mechanics has 
been widely corroborated [7 — 11]. Nonetheless, except 
for the rare possibility to derive analytical solutions F{x) 
of this equation in some cases when V{x) is chosen in a 
very simple and very special form, equation (^ cannot be 
solved otherwise than numerically. Thus, the problem of 
finding a Brownian model to a given quantum-mechanical 
model is relatively difficult.) 

In evolutionary theories of various populations, we may 
use x{t) = Inn(t) — Inn(i), taking n{t) as the number of 
individu als of a certain kind at the time instant t and 
defining Inn(i) as an average value of Inn(i). Since n(t) 
may represent very large numbers, n{t) may be treated 
as a continuous function so that the values of x{t) may 
span the whole set of real numbers. Since mutations are 
random events, a formal stochastic theory of the popu- 
lation genetics can certainly be based on the use of the 
Langevin equation and, correspondingly, of the Fokker- 
Planck equation (|l|). If it is advantageous, we may also 
use the equivalent Schrodinger- type equation (Q). 

Notwithstanding, recently Waxman has shown that 
there is also another mathematical relation between 
the population genetics and the quantum mechanics 
[12]. Waxman's theory concerns a simplified, but well- 
founded, model that we call the Waxman-Peck model. 
(Cf. [13, 14] and references quoted therein.) Wax- 
man's Schrodinger-type equation involves, in the x- 
representation, the 'potential energy' function 



V{x) 



fi exp 



2 2 

(Jrri X 



(6) 



with two mutation parameters, /i > and CTm > 0. The 
'kinetic-energy' operator T in Waxman's equation was 
chosen in the usual form, T ~ — /dx^. 

In the present paper, we will generalize Waxman's the- 
ory. In the momentum representation, Waxman's 'kinetic 



energy' is quadratic, {p\T\p) ^ p^, as it is in the quantum 
mechanics. On the other hand, we assume that {p\T\p), 
being a positive function, need not be quadratic; we only 
require its analyticity along the real p-axis. This vari- 
ability offers further possibilities to model the genetic 
evolution by adequate fitness functions (cf. Section II). 
Obviously, if {p\T\p) > is non-quadratic, the mathe- 
matical relation between the evolution equation of the 
population genetics and the quantum theory becomes 
somewhat more sophisticated than in Waxman's case. 
Namely, when paying heed to the x-representation, we 
have generally to consider a more complicated equation 
than the real Schrodinger equation: in general, our equa- 
tion, with the Hamilton operator E(j)) + V{x), where 
E{p) ^ c(p), is a /unciiona/ differential equation. (Note 
that p — ~ id/dx. The McLaurin development of c(p) 
may involve an infinite number of terms.) For solid state 
theorists, such an equation is familiar as the transformed 
effective mass equation (cf. e.g. the monograph [15] or 
our paper [16]). This equation was invented for envelope 
wave functions of electrons in crystalline solids. (Syn- 
onymously, we may also speak of the real Schrodinger- 
Wannier equation. It is identical with the 'one-particle 
Bloch equation' for the canonical density-matrix with the 
Hamiltonian E{p) + V{x).) 

From the viewpoint of the effectiveness of calculations 
in the present paper, we deem the momentum represen- 
tation better than the ^-representation. Under the as- 
sumption of the smallness of the parameter fj,, we can 
apply the 'plane-wave perturbation theory' (Section III) 
of the density-matrix theory. For a broad class of fitness 
functions, the distribution function of the theory of the 
population genetics can be expressed as a linear expres- 
sion of the mutation parameter ^. 



II. THE WAXMANN-PECK MODEL OF THE 
POPULATION GENETICS 



Let us assume that a large enough habitat (such as 
a given volume of soil) hosts bacteria of a certain kind. 
The habitat yields space, food, moisture, temperature, 
inhibitory substances and other needs for the survival of 
the bacteria in the sense that the total number of the bac- 
teria will never decrease to zero and will never increase 
to infinity. Most of bacteria are free-living microorgan- 
isms multiplying by simple fission. This means that their 
reproduction is asexual. In other words, each bacterial 
individual has only one parent. The typical number of 
bacteria may be huge indeed: one gram of soil may con- 
tain several hundreds million bacteria. Although the bac- 
teria are small organisms — usually 0.3-2 micrometres 
in diameter — their morphology is well distinguishable 
microscopically. 

There are two most frequent shapes of soil bacteria: 
short rods and (slightly deformed) spheres. In both these 
cases, we may characterize each bacterium by its size 
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s. For instance, if the bacterium resembles a rod, we 
define s as the length of the rod. Denoting the average 
of Ins as Ins, we define the phenotypic parameter as 
p = Ins — Ins. Recalling to biology, we consider the 
phenotypic parameter as an inheritable value. If there 
were no mutations in the reproduction of the bacteria, 
all the bacterial individuals in each generation would be 
equally long, i.e. p would be a constant equal to po = 0. 
But then, the mutations — however infrequent they may 
be — account for a dispersal of the value p among the 
individuals, despite the fact that all individuals under 
consideration do still belong to the same biological type. 

Generally, the theory has to respect both mutations 
caused by environmental effects and spontaneous muta- 
tions. For the sake of simplicity, we will consider no 
other than the spontaneous mutations. The spontaneous 
mutations are mainly due to transcription mistakes in 
the replication of the DNA, i.e. in the transmission of 
the genetic information just at the reproduction events. 
Let /i > be the probability of the occurrence of such 
a mutation. In the case of soil bacteria, biologists have 
estimated the values of n between 10~* and 10~^. If a 
mother bacterium is the carrier of the phenotypic value 
p, the daughter bacterium will carry the same value p 
with the probability equal to 1 — /i. (We assume that 
a new-born daughter bacterium grows quickly enough to 
the adult size before becoming mature so that we need 
not distinguish between the size of young and adult bac- 
teria.) If the birth of the daughter bacterium is accom- 
panied with a mutation of the DNA, then there is a non- 
zero probability M{q — p)dq for the possibility that the 
phenotypic value q of the daughter bacterium may lie 
in the interval {q,q + dq), provided that the phenotypic 
value of the mother bacterium was equal to p. Following 
Waxman, we take the function M{p) as a Gaussian, 



M{p) = M{p; am) = 




Here cr^ is the dispersion of values of p. 

Now, to formulate the evolution equation of the pop- 
ulation genetics, we have to introduce the average gen- 
eration time T. Simplifying the problem, we may con- 
sider a discrete time variable as follows. Let the births of 
the bacteria happen at the time instants tn = {n — 1)t, 
n = 1,2,... Then we may say that the bacteria of the 
nth generation live between and i„. Thus, n is the 
generation index. The time discretization is an auxiliary, 
rather formal, mathematical trick which loses its signifi- 
cance if the time t is continualized. For each n, we define 
the distribution function so that <Pn{p)<ip may be 

interpreted as the probability of the occurrence of the 
phenotypic value p in the interval {p, p + dp) in the nth 
generation. The basic problem is to relate the distribu- 
tion function of the generation number n + 1 ('generation 
of daughters') with the distribution function of the gen- 
eration number n ('generation of mothers'). 

Before writing the recurrent formula between the func- 
tions (Pn+i{p) and $n{p), which is our primary goal in 



this section, we have still to mention one important point. 
Even if we have neglected the environmental infiuence 
upon the mutations, we do have to consider environmen- 
tal effects in a Darwinian sense. Namely, we have to re- 
spect that not all bacteria, after their birth, are equally fit 
to survive over the whole generation time r. Only those 
bacteria whose age is equal to r give birth to offspring. 
Some of the bacteria die before becoming mature. These 
bacteria do not take part in producing the individuals of 
the next generation. (However, we assume that even the 
fittest mother bacterium dies soon after giving birth to 
the daughter bacterium. Therefore, we do not include 
the mother bacteria in the number of the bacteria liv- 
ing in the time interval (t„,t„_|_i). The mother bacteria 
have been included in the number of the bacteria living 
in the time interval {tn-i,tn)-) The fitness of the bacte- 
ria to live in their environment until their maturity can 
be modeled by a non- negative function w{p). Requiring 
that 

< w{p) < 1 , (8) 

we may give the function w{p) a probabilistic meaning. 
We assume that a new-born carrier of the phenotypic 
value p has the chance to live until the maturity with 
the probability w{p). The number of mature carriers 
of the phenotypic value p from the interval (p, p + dp) 
in the nth generation is proportional to w{p)$n{p)- To 
determine the shape of the function w{p) should be a 
matter of thorough biological investigations from case to 
case. We suppose, as Waxman and Peck did, that w{p) 
behaves analytically around the value po = and that 
this value corresponds to the maximum value of w{p). 
(Apparently, the fittest bacterial individuals are those 
whose phenotypic parameter p is equal to the average 
value p. However, p = 0.) 

Waxman and Peck have chosen the function w{p) in 
the special form 

w{p) = w{0) exp(- 7p^) , < w{0) < 1 , (8a) 

assuming that < 7 <C 1. (In fact, expression (8a) 
defines the Waxman-Peck model. The value of ii;(0) is 
insignificant since the distribution functions $n{p) are 
independent of w{0).) 

There are, of course, many other possibilities to model 
w{p) by slowly varying functions with the maximum at 
Po = 0. These functions need not tend to zero if \p\ ^ 00. 
(The fitness function has been defined as a probability, 
not as a probability density!) In order to illustrate how 
the theory may depend on the choice of the function w{p), 
we will treat, in addition to the Waxman-Peck model, 
also an alternative model. Our model (considered as an 
example) is defined by the fitness function 

w{p) = w{0) {1 - 7[1 - exp(- ap^)]} . (8b) 

Here we suppose that < 7 <C 1, admitting that a > 
need not be a small number. Expression (8b) tends to 
w(0)(l - 7) > if IpI ^ 00. 
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The distribution function $n+iip) of the generation 
of daughters is determined by two contributions from 
the generation of mothers. The first stems from the 
births without mutations. The phcnotypic value p is un- 
changed at such births and the probabihty of the occur- 
rence of such births is equal to 1 — /i. The first contri- 
bution to ^?„+i(p) is proportional to (1 — fi)w{p)(l>n{p). 
The second contribution to <?„+i(p) is proportional to 
fj, J^^dq M{p — q;crm)i'j{q)^n.{q)- The interpretation 
of this expression is clear: if the birth of the carrier 
of the phenotypic value p is accompanied with a mu- 
tation, we have to consider mature individuals, allow- 
ing all possible phenotypic values q of potential mothers. 
To exhaust all such possibilities, we have to integrate 
M(p— q]ain)w{q)<l>niq) with respect to q. Since both 
$n(,p) and (!>n+i{p) are probability densities, we have to 
require that 



dp <Pn{p) 



dp <Pn+i{p) = 1 



(9) 



Therefore, we write the equality 



(1 - n)w{p)$n{p) + M{p - q- am)w{q)^n{q) 



(10) 



The denominator in the r.h. side of equation (|10|) war- 
rants the fulfillment of condition (||). Since 7 is small, 
it is convenient to introduce the complementary function 

c(ji) to w{p)/w{0): 



c{p) = 1 



w{p) 



(11) 



In the case of the Waxman-Peck model, 



c(p) = 1 - exp(- 7p2) 



(11a) 



whilst in the case of the model defined by function (8b) , 



c(jp) 7 [1 - cxp(- ap^)] 



(lib) 



From the viewpoint of biology, the smallness of 7 implies 
that the comparison of the survival fitness of the major- 
ity of the bacteria with the survival fitness of the fittest 
bacteria should not reveal too conspicuous differences. 
When using the function c(p), we can rewrite formula 
(0) in the for m 



^n+l{p) 



(1 - m) [1 - c{p)] jp) + ^l f^^ dq Mjp - q; am) [1 - c(g)] ^„ (g) 
l-/r^dgc(g)f„(g) 

I 



(10') 



With realistic values of p around po = 0, the values of 
7P^ are small. Thus, in the case of the Waxman-Peck 
model, the values of c{p) are also small and 



C{p) = 7p2 + 0(^2) 



(ll'a) 



On the other hand, in the case of the model defined by 
function (8b), we have to keep expression (lib) intact 
since a need not be a small parameter. 



We may take advantage of the possibility to neglect all 
terms of the order of magnitude of 7^, as well as of 7/x. 
So we write 



l-/!ldg c{q)Mq) 
and 



= 1 



dqc{q)$n{q) 



J 



Pn+l[p) = Mp) 



c{p) 



dq c{q)<Pn{q) 



I 



/OO 
dgM(p- g;crm)<l>„(g) 
-00 



(10") 



Now, in the approximation neglecting the terms sym- 
bolized by the dots, we are ready to go over into the 
formalism employing the continual time variable t, real- 



izing that the value of the generation index n may be 
high. Typically, the generation time r of soil bacteria is 
about 20 minutes. This means that after elapsing hun- 
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dred days, the genetic information passes over more than 
seven thousand generations of the bacteria. If n ^ 1, 
we may identify $n{p) with <l>{p,t) and approximate the 
difference <?„+i(p) — ^n(p) as the time derivative: 

<Pn+i{p)-Mp)=r^^^ + ... (12) 



Thus, we can rewrite equation (10") in the approximate 
integro-differential form 



at 



c{p)- / dqc{q)<P{q,t) 



r 



/oo 
dqM{p^q;<Jm)'P{q,t) 
-OO 



(13) 



This equation was derived in [12] (where, however, c(p) 
was approximated as 7P^). Evidently, equation (|l^) is 
nonhnear. Fortunately, this nonlinearity does not mean 
a serious problem since we may employ the substitution 



Hp,t) 



ip{p, t) 



and require the validity of the equation 



(14) 



^^^^ = --c(p)<^(p,i) + - / dqM{p~q;am)^{q,t)- 
at T T 

_ (15) 
Equation (|l5|) is linear. After integrating it with re- 
spect to p, we obtain the equation d/dt dq (p{q, t) = 
^ ( V't) dpc(p)¥3(p, t) + {^i/ t) dqip{q, t) and when 
substituting expression ( [T^ ) for (p{p,t), we arrive at the 
identity 



d_ 
di 



dqip{q,t) 



dqipiq,t) 



(16) 



dpcip)^{p,t) + ■ 

The differentiation of expression (|lj) gives the identity 
dip{p,t) _ 



dt 



dHP,t) 
dt 



f°° d r°° 

J dg^(g,t) + <l>(p,i) — J dqip{q, 



t). (17) 



When equalizing the r.h. sides of equations ( |15| ) and (|lj 
and when respecting identity (|l^), we obtain equation 
(13) for the function <P{p,t). Thus, instead of directly 
solving equation ([T^), we may solve Waxman's equation 
( p^ ) at first. This task, as we will show in Section III, is 
not difficult. If the function <l>{p, t) obeys linear bound- 
ary conditions, the function ip{p, t) has to obey the same 
boundary conditions. Wc will simply assume that 



$ip,t) and ip{p,t) if \p\ 



(18) 



It remains still to discuss the initial condition. Whichever 
initial function 



(19) 



${p,0) = <Po{p) 



is chosen, the solution <P{p, t) for t > of equation (|lJ 
is unique. Since equation (Qq) is linear, we may mul- 
tiply ip(p,t) by an arbitrary constant A. If ip{p^t) gives 
the function <^(p, t) then A(p{p, t) does also give the same 
function 4>(p,t). Therefore, we may choose the integral 
dq(p{q, 0) (which is a constant) equal to unity. Then 
formula ( p^ ) and equality (|l^) give us the initial condi- 
tion 



ip{p,0) = Mp) 



(20) 



for the function ip{p,t). 

If <?o (p) is an even function, equation ( [l5| ) implies that 
the function ip{p,t) is also even in the variable p and 
<I>{—p,t) — ${p,t) at all times t > 0. In this case, the 
mean value of p is an invariant in time (i.e. a constant): 



/oo 
dp p <l>{p,t) = . 
-oo 



(21) 



Equation (|T^) is formally the same as the Schrodinger- 
Wannier equation in the momentum representation. It 
can easily be Fourier-transformed. We define the func- 
tion 



ipix^t) 



1 



/2Tr 



dp e:>q){ipx)(p{p,t) . (22) 



This function is the solution of the functional differential 
equation 



dip{x, t) 



1 



exp 



^2 ^2 



ilj{x,t) . 



T \ 2 

In Waxman's approximation, equation (^3|) reads: 
dtpix,t) 7 d^ipix^t) 



(23) 



dt 



dx"^ 



6 



cxp 



v2 ^2 



(23a) 



If c{p) is taken in the form of expression (lib), the func- 
tional differential equation for 'il}{x,t) reads 



dij}{x,t) 7 



exp 



° dx'- 



■ipix,t) 



v2 ^2 



ip{x,t) 



(23b) 



III. THE PLANE- WAVE PERTURBATION 
THEORY 



Instead of solving equation (23a) or equation (23b) and 
carrying out the integration 



vip, t) 



1 



%/2^ J- 



dx exp(— ipx)'il){x, t) , (24) 



we prefer to calculate the function Lp{p, t) directly. Defin- 
ing the 'potential-energy operator' V{p) 

V{p)^{p,t)^- f dqMip-q;amMq,t) , (25) 

let us write equation (|l^) in the form 

d(p{p, t) c{p) 



dt 



ip{p,t) + fiVipMp.t) (26) 



and define the Green function G{p, t;po) of this equation. 
Employing the Green function, we write f(j>, t) (for t > 
0) as the integral 



dpoG(p,t;poMpo,0). (27) 



The initial function ip{p, 0) has been defined by equality 
(pO|). The Green function itself obeys the equation 



dG{p,t;pa) ^ c(p) 
dt T 



Gip,t■,po) + ^lVip)G{p,t;po). (28) 



According to equality ([27|), G(p, i;po) satifies the initial 
condition 



G(p,0;po) S{p-po) 



(29) 



In the special case when c{p) is approximated by the 
quadratic function, equation (^8|) is formally identical 
with the Bloch equation for the one-particle canonical 
density matrix Cp{p,pa) in the thermodynamic equilib- 
rium. In the case of a general function c(p) we have to 
speak of quasiparticles with a non-parabolic dispersion 
law. The function C-i-o(PiPo) is equal to 5{p — po) for 



quantum-mechanical reasons. When transforming equa- 
tion ( p8|) into the a;-representation form, one observes 
that the 'potential energy' corresponds to a well: it is an 
inverted Gaussian (cf. expression (^). 
We can derive G{p, t;po) as the series 

oo 

G(p,t;po) = 5Z^'^j(P'^;Po)- (30) 

The zero-order term is the solution of the equation 

dGo{p,t;po) c(p) 

^7 = Go[p,t;poj (31) 

Ot T 



with respect to the condition 

Go{p,0;pq) = 5{p-po) 



(32) 



When solving equation (|31|), we obtain, for t > 0, the 
function 



Go{p,t;po) = (5(p-po)exp 



c{p)t 



(33) 



With this function, we can write down the integral form 
of equation (p8|): 

G{p,t;po) = Go{p,t;pQ) 
+ IJ. dti dpi G{p,t ~ ti;pi)V(pi)Go{pi,ti;pQ) . 

Jo J-oo 

(34) 

This is a Dyson-type series (cf. e.g. [17].) The first-order 
term in this series (linear in /x) reads: 

fiGi{p,t;po) = 

nt POO 

dti dpiGo{p,t-ti;pi)V{pi)Goipi,ti;po) 

Jq J -oo 



pt /"OO 

■ / dti dpi Go{p,t - ti;pi) 

Jo J -oo 



X / dqM{pi- q]am)Go{q,ti;po) ■ (35) 
<y —oo 

After inserting expressions (0) and ( ^ ) here, we obtain 
the function 



nGi{p,t;po) = 



fif 1 



T \27rCT- 



1/2 



exp 



{P - Pof 
2^2 



[\u f c{p){t ~ ti) + c{po)ti 
X / d^i exp ' 



After performing the integration with respect ti, we ar- 
rive, respecting formula (M), at the final result 
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.J, N exp[- c{po)t/T] ~ exp[- c{p)t/T] 

c{p) - c{pq) 

In this same way, we could also calculate higher-order 
terms (i.e. the terms proportional to ^l^ with j > 1) 
in series (^). We expect, however, that higher-order 
terms are negligible, since the mutation probability is, 
as biologists have proved in their extensive studies, very 
small. 



at all times t > Q. Hence, in the linear approximation 
with respect to /i, we have got the function 

^(p, t) ^ 5{p) + mM(p; am) ^ ~ '"P^^'^)'^^^'^"^ • (41) 

The only problem that we have still left unsolved is the 
calculation of the integral 

/_^dp^(p,t) = l+,(^) ' X 



IV. DEVELOPMENT OF THE PHENOTYPIC 
DIVERSITY IN A POPULATION WHOSE 
INDIVIDUALS ARE INITIALLY EQUAL 



If all individuals of a population are initially, at the 
time = 0, carriers of the same phenotypic value poj 
the initial distribution function <?o (p) is equal to the delta 
function: 



Mv)^5{p) 



(37) 



(po = p ~ Q according to our definition of the phenotypic 
parameter p.) In regard to identity (|20|), equation (|2^ ) 
allows us to assert that 



(^(p,t) = G(p,i;0) 



(38) 



Since c(0) = (cf. expression (pi])), formulae (|3^) and 
(pq) imply, respectively, that 



Go(p,i;0)=,5(p) 



and 



/iGi(p, t. 0) = ^iM{p] cTm) 



1 — exp[— c{p)t/T] 
c{p) 



(39) 



(40) 



dp exp 



p^ \ I — exp[— c{p)t/T] 



2a^ 



z{p) 



Recall that, according to (|l4|), 

Hp,t) = N{t)^{p,t) 



where 



N{t) 



/OO 
dpip{p,t) 
-OO 



(42) 

(43) 
(44) 



We will calculate the function N{t) approximately, as- 
suming that < (Tm < 1. (In fact, it is probable that 
Cm ^ !•) The most relevant values of p contributing to 
the value of the integral in the r.h. side of formula (E2) 
lie in the interval (— Urn, (Jm)- 



A. Distribution function <l>(p,t\Q) in the model 
where c(p) = 1 — exp(— 7p^) (the Waxman-Peck model) 



Since < 7cr^ <C 1, we may use, when calculating 
integral (p2[), the approximation expressed by formula 
(ll'a). Thus, 



dp(p{p,t) w 1 + ^ 



1/2 



dp exp 



p"^ \ 1 — exp(— 7p^i/T) 



2(7^ 



M f 1 



T V 27r(Ti 



1/2 /.CXD 

2 , / dii / dp 

m/ -/o i-oo 



exp 



2fT2 



7*1 \ 2 



r 



After carrying out the integration with respect to p, we 
obtain the simple result 



27cr^^t 



2 +\ 1/2 



/ , M /■* dii 

dp<y3 p,i W 1 -I- - / — — 

T Jo (l + 27criVi/r)i/^ 



Hence, according to formula (|TJ), we obtain the distri- 
bution function 
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#(p,t;0) 
Nit) 



d(p) + -M{p;(7m) 5 

7 P 



(43a) 



with the normalizing coefficient 



N{t) 



70-m 



1/2 



(44a) 



From the probabilistic viewpoint, the function N{t) is 
well understandable. When counting all bacteria living 
at the time instant t, we have to distinguish whether 
they are carriers of the original phenotypic value po — 
or whether they carry other values, p ^ 0. Since 
Jle ^(P' ^) ~ ^(^) i^^ ^ ~^ +0), we may say that a ran- 
domly chosen bacterium may be the carrier of the value 
Po = with the probability equal to N{t). U t 00, 
the probability N(t) decreases towards zero. However, 
this decreasing — the process influenced both by the mu- 
tations and by the fitness of the bacteria to live in their 
environment — is slow. Indeed, let us take 7 = 0.02, 
CTm = 0.05 and ^ = 10"^. Then fi/{j(T^) = 0.2 and 
2ja'fjyt/T — 1 for the generation number t/r — 10"'. If 
these values of 7 and am, together with the value 20 
minutes for the generation time r, may be taken as re- 
alistic for some soil bacteria, the total time t comprising 
the lifetime of ten thousand generations of these bacteria 
equals about five months. If the time t is roughly ten 
times (or more than ten times) shorter, formula (44a) 
can be simplified: 



we may approximate 1 — exp(— ap^) by unity. Then we 
obtain the simple result 



/ dpipip,t) « 1 + ^ 
7 

Correspondingly, if 017^-^ 3> 1, then 



1 ( 

1 — exp I — — 



N{t) { Siz) + ^ M{p; am) 



where 
Nit) 



l-exp( - ^ 



exp 



T 



(43b) 



(44b) 



In the short-time approximation, formulae ( |45| ) and (45') 
are equally valid as in the case A. Note that expression 
(43b) for the distribution function (l>{p, t) would be cor- 
rect if c{p) = 1 — w{p)/w{0) might be approximated by a 
small constant 7 > 0. In this case, N{t) may again be ap- 
proximated as 1 — fit/ J at short enough times. However, 
at 00, then N{t) does not tend to zero (in contrast 
to the case analyzed in the preceding subsection): 

lim N{t) = — . 
t ^ 00 7 + M 



V. CONCLUDING REMARKS 



N{t) 



^1 if 



27<i 



< 1 



(45) 



As a rule, the mutation probability fi is smaller than 
2jam- Thus, we may write 



N{t) 



^ if 



27c^mi 



< 1 



(45') 



Distribution function #(p, t;0) in the model 
where c(p) = 7[1 — exp(— ap^)] 



Now we consider a small parameter 7 (0 < 7 <C 1) 
and another parameter, a > 0, which need not be small. 
Only if a <C l/f^, we may accept the approximation 
c(p) K, ^ap? and there is no essential difference from the 
Waxman-Peck model, only 7 is replaced by 7a. 

Otherwise, if atr^ is comparable with unity, the inte- 
gration of the function (/^(p, t) with respect to p is much 
more complicated but can be accomplished explicitly. (It 
is presented in Appendix.) 

Here we confine ourselves to discussing what comes 
about if acr^ 3> 1. Essentially, under this condition, 



In the present paper, we have focused attention on 
the importance of the fitness function w(j)) in the theory 
of the population genetics. Assuming that < c(p) = 
1 — iy(p)/w(0) <C 1, we have essentially followed Waxman 
and Peck who derived the distribution function #(p, i) 
of the population genetics as a functional of a function 
Lp{p,t) (cf. expression ([ij)) satisfying a linear integro- 
differential equation (cf. equation (^5|)). However, in 
contrast with paper [12] where c{p) was approximated as 
7P^ with some small parameter 7 > 0, we emphasize that 
c{p) may be chosen from a wider class of functions. In 
particular, we have dealt with the model defined by the 
function c{p) ~ ^[1 — exp(— ap^)]. 

We have calculated the distribution function as a series 
with respect to the mutation probability fi. Our iteration 
scheme for calculating the Green function G{p,t;pQ) of 
the equation for (p{p, t) has been used in the same manner 
as in the density-matrix theory. 

The replacement of c{p) by E{p), t hy (3 (with h — 
1) and Go{p,t;po) by the unperturbed canonical density 

matrix Cp*\p,po) yields the equation 



9Cf 



iP,Po) 



d(3 



E{p)cf{p,po). (46) 
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With adequately chosen function E{p), this equation 
may concern conduction electrons in a homogeneous non- 
degenerate semiconductor. (Since E(p) is not equal to 
the kinetic energy of an electron in vacuum, we may in- 
terpret the conduction electrons as quasiparticles defined 
by the dispersion law E — E{p).) 

Our second remark concerns the analogy with the dif- 
fusion theory. The Fourier transform of the function 
G{p,t;po) (multiplied by a constant) can be interpreted 
as the concentration C{x,t;xQ) of diffusants which all 
were initially, at the time to — 0, localized in the point 
xq. In the approximation of the present paper, we may 
generally write the equation 



dCa{x,t;xa) 
dt 



1 



d 



The Poissonian process consists of equal delta-pulses: 
all the pulses are taken with the same amplitude a. 
The average frequency of these pulses is equal to 7/r. 
Clearly, we consider a multiplicative stochastic process 
x{u) (0 < u < t). If Dq — 0, then the probability 
density Po{x^ t; xq) = {6{x — x{t))) is the fundamental so- 
lution of equation ( 4S ) . Alternatively (as we have shown 
in [18]), equation (4S) can be written in the equivalent 
integro-differential form: 



dPo{x,t;xo) _ 7 
dt ~ T 



dx' 



1 



(27ra)i/^ 



exp 



{x - x'f 



2a 



6{x — x) 



Po{x,t;xo) + Do 



d'^Po{x,t; Xq) 



dx'^ 



(49) 



exp 



(Jyyi X 



C{x,t; xo) 



(47) 



If c{p) — 7p^, the concentration C{x,t;xQ) obeys the 
usual diffusion equation with the diffusion coefficient D = 
7/t. If c(p) 7^ 7p^, the diffusion is anomalous. In 
any case, the positiveness of the 'potential-energy term' 
means that equation (^) involves a creation of diffusants. 

If /i = 0, we observe that No = dx Co {x, t; xo) is a 
quantity not varying in time. Therefore, we may define 
the probability density Po{x,t;xo) = Co{x,t;xo)/No and 
put the theory on an equal footing with the theory of the 
Brownian motion. 

If c(p) = 7p^ and fi — 0, we may write down the 
Langevin equation x{u) = (27/r)^/^/(it) for the stochas- 
tic paths x{u) {0 < u < t) which all start from the 
common point x{0) = xo at the time instant uo = 0. The 
value of the end-point a;(t) = a; at a given time instant 
t > may be arbitrary and Po{x,t; xo) = {S{x — x{t))). 
In the terminology of the theory of stochastic processes, 
x{u) is the Wiener process. 

But then a natural question arises: which stochastic 
process does correspond to the case when the fitness func- 
tion w{p) is modeled by function (8b) with which we have 
exemplified our problem? About ten years ago, we dealt 
with the equation 



dPo{x,t; Xo) _ 7 
dt ~ T 



exp a 



1 



Po{x, t; Xo) 



In the case when Do = 0, equation (^) was employed 
by Laskin [20] in a theory of the channeling of high- 
energy particles in crystals. (The channeling occurs when 
a ray of equi-energy particles bombarding a crystal is col- 
limated very precisely in a favorable direction.) 

In the framework of the diffusion theory, we may con- 
clude that the parameter 7 of the theory of the popula- 
tion genetics corresponds to an environmental noise. If 
7 = 0, the noise is absent. 

Section IV of the present paper has been devoted to 
the problem of the evolution of a population in which 
all individuals are initially equal, being the carriers of 
the phenotypic value po — 0- The distribution function 
<?(p, t; 0) of the population is a sum of a sharp delta- 
function component, N{t)S{p), and a blurred component. 
Similarly as in the thermodynamics, we may distinguish 
two phases in the population at any time t > 0. Let 
us denote them as phase S and phase B. The phase S 
consist of the carriers of the initial phenotypic value po = 
0. The phase B consists of the individuals carrying the 
phenotypic values p ^ po- In the Waxman-Peck model 
(cf. expressions (43a) and (44a)), the probability N{t) 
tends to zero if t —^ 00. Therefore, we may say that the 
phase S dissolves gradually in the phase B. In the model 
with the fitness function w{p) defined by expression (8b) 
(or by another similar expression), the probability N{t) 
does not tend asymptotically to zero: this model predicts 
that both the phases S and B may coexist if t — > 00. 



(Cf. equation (45) in [18]; see also [19].) Equation 
( ^ ) corresponds to a stochastic process with paths x{u) 
defined by the stochastic equation x{u) = [{2Do)^^^ + 
aJ2jH'^ ~ %)] fi'^)^ where f{u) is the standard zero- 
centered Gaussian white-noise function and where the 
sum represents a point process in which uj are ran- 
dom time instants distributed in the Poissonian way. 
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APPENDIX 



Thus we have obtained the distribution function 



In the model where c(p) = 7 [1 — exp(— ap^)], we have 
to manage the function 

1 — exp[— c(p)f/r] 
c{p) 

1 — exp { — 7[1 — exp(— ap^)]t/T} 
7[1 — exp(— ap'^)] 



#(p,i;0) = iV(i) 



Sip) 



H 1 - exp { - 7[1 - exp(- ap'^jt/r ' 

+ - M{p; am) z -, Trr\ 

7 1 — exp(— ap-'t/T) 



where 



if 

T- Jo 



dti exp 



7[1 - exp(- ap^)] h 



dti exp ( — — exp — exp(— ap ) ] . 
T Jo \ r J \ T ) 

Wc have to calculate the integral 
I" dv^{p,t) = 1 + ^ ^*dti exp ^) /(ii) , 



where 



1/2 ^00 



X exp ( — exp(— ap ) 



When developing the second exponential in the MacLau- 
rin series, we obtain the following sum of the Laplace 
integrals: 



1/2 00 



dp exp 



1 



2(7, 



m 



Hence, 



^.^ j\T^l + 2jaa^)y^ ■ 
In this way we have obtained the result 



E 



dti exp 



7t] 



•-1 • 



The integral in the r.h. side of this equality is easily 
calculable: 



/c=0 



^ -.fc^fe 



E 



^^7^"(l + 2iaa^)V2 



1 — exp 



7t\ 7'^t'^ 



fe=o 
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